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We present a stochastic approach for ion transport at the mesoscopic level. The description takes 
into account the self-consistent electric field generated by the fixed and mobile charges as well as the 
discrete nature of these latter. As an application we study the noise in the ion transport process, 
including the effect of the displacement current generated by the fluctuating electric field. The 
fluctuation theorem is shown to hold for the electric current with and without the displacement 
current. 

I. INTRODUCTION 

The study of electrolyte solutions finds its origin in the seminal works of Nernst [TJ dj and Planck [31 0]. Ion 
transport in electrolytes is described as arising from a diffusive part due to the concentration gradients and a drift 
part due to the external electric held imposed to the system. A more accurate description is obtained if the electrical 
held is not only generated by external means but also self-consistently incorporates the contributions of the local 
deviations from electroneutrality [SJ [5] . Under the assumption that the electric held propagates instantaneously - 
the quasi-static limit of the Maxwell equations - this problem is known as the Nernst-Planck-Poisson (NPP) problem 
H] . Here, the electric held self-consistently arises from the distribution of charges as described by the Poisson 
equation. In this case, the transport process may present an additional contribution due to the temporal variations of 
the electric held, known as the displacement current [9THT] . This approach has also played a basic role in the theory of 
other systems: Indeed, the NPP equations, combined with statistical arguments, lead for instance to the description 
of semiconductor p-n junctions [T2l IT5] . 

On the other hand, ionic solutions exist in many natural or artificial systems. Many biological processes crucially 
depend on the transport of ions, e.g., between intra- and extracellular solutions [T3KH]. Besides, solid-state nanopores 
are fabricated for the study of ionic current fluctuations [TB]. In these systems, the transport of ion takes place in 
different geometries, either with the complexity of channel proteins or more regular as straight channels in the case 
of nanopores. Biological ion channels [17], as well as physico-chemical systems such as nanopores [18], nanofluidic 
diodes [IS], or nanostructures studied by impedance spectroscopy [20]) are often described by the NPP theory. Yet, 
they are of mesoscopic sizes and, at these scales, the motion of ions in solution is subjected to molecular fluctuations 
and presents a stochastic behavior. These random fluctuations can be successfully described at the mesoscopic level 
as Markovian stochastic processes [2"D - f2"5] , Non-Markovian discrete-state models have also been considered to take 
into account the memory effects resulting, in particular, from the lumping of diffusion [2"B] . The description to be 
adopted depends on the geometry. In short heterogeneous channels, transport tends to proceed by jumps between 
discrete states corresponding to wells in the free-energy potential. In long homogeneous channels such as cylindrical 
nanopores, diffusion is not interrupted by barriers except at entrance and exit. In any case, the master equation ruling 
these stochastic processes should reduce to the evolution equation for the charges densities at the macroscopic level. 
In this respect, one of the fundamental problems is to incorporate the long-range Coulomb interaction between the 
electric charges in the description. 

In the present work, our purpose is to present a description of ion transport which is consistent with the laws of 
both electricity and statistical thermodynamics. The key point is that the long-range Coulomb interaction deeply 
influences the fluctuations of the particle and total currents. The proposed model describes the spatial distribution 
of the discrete number of ions in a long channel. The ions undergo random jumps due to the thermal agitation in a 
self-consistently generated electrical potential landscape. This description takes into account the stochastic aspects 
of the time evolution as well as the self-consistent electric field generated by the charge distribution. In this regard, 
it offers a computationally favorable alternative to more expensive techniques such as Brownian dynamics [271 [35] . 
This approach allows us to study various aspects of ion transport, revealing, e.g., the long-ranged spatial correlations 
of the potential fluctuations inside the channel. We further study the fluctuations present in the transport process, 
including the contributions of the displacement current. Also, we show that the fluctuation theorem [2"9"H5fj] describing 
the large-deviation properties of the current holds in this case as well. 

The paper is organized as follows. In Sec. [TTJ we introduce the macroscopic description of ion transport in terms 
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III The fluctuations in the 



of the Nernst-Planck-Poisson equations. The stochastic description is presented in Sec. 
ion transport process are studied in Sec. IV In Sec. |VJ the symmetry of the fluctuation theorem is shown to hold for 
both the particle and electric currents. Conclusions are drawn in Sec. [yT] 



II. MACROSCOPIC EQUATIONS 

We consider a one-dimensional conduction channel extending from x — to x — I. An important class of electrolytic 
systems are those which contain a homogeneous distribution of immobile ions [37] . We therefore consider a model 
with a single mobile species of charge e and density n being transported in a channel presenting a fixed ion density 
?i_ of opposite charge — e. The charge density inside the system is thus given by 

p = e(n — n_) . (1) 

The case of several fixed and mobile species can be treated in a similar way. Expressing the electric field £ in terms 
of the electric potential <I> as £ = — V$, the problem is ruled by the coupled diffusion and Poisson equations: 

d t n = V • (DVn + ^enV$) , (2a) 

V 2 $ = --(n-n_). (2b) 
e 

D is the diffusion coefficient of the mobile species and [i is the mobility coefficient which is given by Einstein's relation 
H = D/k^T. The dielectric fluid is here homogeneous throughout the channel and has a dielectric constant of e. 



Equation (2a I is the conservation equation for the particle density, dtn = — V • J, where J is the particle current 



density. The Poisson equation (2b) in turn gives the electric potential as a function of the charge density. This 
coupled system ([2| forms the Nernst-Planck-Poisson equations. They correspond to the quasi-static limit of the 
complete Maxwell equations. On the other hand, the total current density reads [HI E] 

I = eJ + e-, (3) 

the last term being the displacement current density arising from temporal variations of the electrical field. The 
displacement current vanishes in a stationary state but contributes, for example, when the system relaxes toward 



the stationary state. Note that, by virtue of the Poisson equation (2b), the total current density is divergence free, 
V • I = 0, at all times. We further notice that the experimentally measurable quantity is given by the total current 
[ED- 

The channel is in contact with two reservoirs of particles maintained at fixed concentrations and electric potentials. 
Equations ^ are thus supplemented by the boundary conditions 

n(0) = n L , n{t) = n R , (4a) 
$(0) = $ L , $(£) = $ R , (4b) 

determining the concentrations on the left (L, x = 0) and right (R, x = £) boundaries as well as the external potential 
difference V — $l — $r applied to the channel. 

For time-independent boundary conditions, the system evolves toward a stationary state where d t n = so that 
the particle current density J is constant in time and in space. For the one-dimensional channel here considered, the 
stationary density satisfying the boundary conditions Q is expressed as 

n(x) = n L e^'^ - - e*< x > / e'^ dy , (5) 

along with the particle current 

T n n R e"^ - n L e-^ 

3 = ~ D Vi T7~~ . (6) 

where we introduced the dimensionless potential <j> = —e^/k^T to simplify the expressions. This potential obeys the 
dimensionless Poisson equation 
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where 

y Th- 
is the Debye screening length of the system [35] . 

The stationary state is a state of thermodynamic equilibrium if the particle current vanishes, i.e., when we have 

(0L-0R.)c q = ln^ or (<& L - * R ) oq = ^ln^. (9) 
n R e n L 

This voltage difference is known as the Nernst potential [H[2]. If this condition is not satisfied we are in a nonequi- 
librium situation characterized by the presence of a non- vanishing ionic current given by Eq. ([6]). In the stationary 
state, the associated total current reads 

I = eJ (10) 

since the displacement current vanishes in this case. However, at the mesoscopic level, fluctuations in the charge 
distribution generate a fluctuating displacement current, as studied in Sec. |IV| 



III. STOCHASTIC DESCRIPTION 

The study of noise and molecular fluctuations at the mesoscopic scale is successfully accomplished in terms of 
a Markovian random processes. Indeed, the master equation is known to describe the fluctuations down to the 
mesoscopic scale [2"lTf2"3"] . 

We here introduce a stochastic model for the distribution of ions in the channel that incorporates the self-consistent 
electric field generated by the fluctuating distribution of ions. The channel is divided into L cells of volume AV, cross 
section a, and length Ax = l/L centered at the positions Xj = (i — 1/2)/ Ax (i = 1, . . . , L). Each cell contains a 
discrete number of mobile ions Ni and a constant number N_ — n_AV of fixed ions. The cells and L + l correspond 
to the external reservoirs maintained at fixed concentrations so that their particle numbers Nq and Nl + \ remain 
constant in time. The dimensionless electric potential fa — — e^/fceT is defined on each cell as well. It obeys the 
Poisson equation ^ which is discretized according to 

fa+i - 2 fa +fa-i If Ni 



Ax 2 A 2 ViV_ 



1 (11) 



with the boundary conditions fa = fa, and 4>l+i = fa\,. Equations (11) form a linear system that must be solved 
at each time the particle distribution changes. Its exact solution is given in Appendix [XJ The extension of the 
model to several (positive or negative) ion species and inhomogeneous distribution of fixed ions and permittivity is 
straightforward. 

The state of the system is determined by the number of ions Ni in each cell. The evolution equation for the 
probability distribution P(Ni, Nl) to observe a configuration {Ni} of particles is ruled by the master equation 

d L 

jP{Nx, N L ) = Y J [ W +i (...,Ni + 1, N i+1 - 1, ...)P(..., Nt + 1, N t+1 - 1, ...) 

i=0 

- W + i(...,Ni, N l+1 ,...)P{..., N u N i+1 , ...) 

+ W-i(..., Ni - 1, N t+1 + 1, ...)P(...,Ni - 1, N l+1 + 1, ...) 

- W-i{..., Ni,N i+1 , ...)P(...,Ni, N i+1 , ...) ] (12) 

where W±%{-) denotes the transition rate between two configurations of the system. The transition ±i changes the 
configuration from (..., Ni, N i+ i, ...) to (..., iVj =p 1, A i+1 ± 1, ...). These transitions rates are supplemented with the 
boundary conditions that No and Nl + \ take fixed values. 
The transition rates can be expressed as 

W +i (...,Ni,N i+1 ,...) = V(AI7 ili+ i)JVi (13) 
W-i{...,Ni,N i+1 ,...) = i>(AU i+lti )N i+1 (14) 
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with the function 



where D is the diffusion coefficient and /3 = (fee? 1 ) -1 the inverse temperature. The prefactor D/Ax 2 allows us to 



recover the appropriate evolution equation (2a) in the macroscopic limit, L — > oo 39J . The transition rates are 
proportional to the number of mobile ions in the cells and depend on the free electrostatic energy difference AC/^+i 
as a result of the transition event: 

AU iti+1 = ~(Vi + V!). (16) 



In Eq. ( 16 ) , Vi and V- are the voltage drops across the cells i and i + 1 before and after a transition event, respectively. 
This expression corresponds to the change in electrostatic energy resulting from the transition event, as shown in 
Appendix [A] Consequently, Eq. (16) can be expressed as 



AUi, i+1 = |($ i+1 - + - (17) 



in terms of the electric potential calculated for the current configuration, $j = <f>j(Ni, Nl), and the electric 

5d: 

$; = $;•(..., NiTl,N i+l ± !,...)■ (18) 



potential that would occur if the corresponding transition ±i were performed 



The potential $>j can be expressed as 

d^^efcr 1 ) iefc- 1 ) (19) 

\ / 7,2 \ / 7,2 + 1 



in terms of the matrix C obtained from the solution of the discretized Poisson equation (see Appendix | A 
The function ip satisfies the identity 

ip{AU) = i;(-AU)e-P AU (20) 



which guarantees that detailed balance is fulfilled at equilibrium [40-42]. On the other hand, the nonequilibrium 
constraints or affinities driving the system out of equilibrium [43] can be identified in the stochastic description 
thanks to a construction put forward by Schnakenberg [22] . according to which the affinities are obtained from the 
cyclic trajectories in the forward and backward directions. In our case, as verified in Appendix [B] the macroscopic 
affinity is readily identified by considering trajectories involving the transfer of an ion from one reservoir to the other, 
yielding 



A = In -^^e A * = In ^e A * . (21) 

Here, A<p = 0r — 0l and we note that this affinity only involves macroscopic quantities as it should. At equilibrium 
the affinity vanishes, ^4 cq = 0, and we recover the macroscopic equilibrium condition 

The stochastic process is simulated as follows. Random trajectories of the system are obtained using Gillespie's 



algorithm [44], which is known to reproduce the statistical properties of the master equation (12). At each random 



jump, we have to recalculate the electric potential in the channel by solving the system (11). The simulations are 
performed with the following parameters. Each cell contains a number iV_ = 25 of fixed ions. By an appropriate 
choice of the time unit, we can impose for example D/Ax 2 = 1. Using the dimensionless electric potential <j), the 
only remaining dimensional parameter is the Debye length A. We form the dimensionless quantity £/X, i.e., the ratio 
between the channel and Debye's length, for which we choose the value lj\ = 50. For a Debye length of 7nm this 
would correspond to a channel length of 350 nm. 

The link with the macroscopic description is established in the limit where the volume of the cells AV vanishes in 
which case the concentrations are recovered as n(xi) — (iVj) /AV. In Fig. [l^,, we depict the average number of ions 
along the channel, estimated by the time average 

(A, ; ) = lim i / Ni(t)dt (22) 
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FIG. 1: (a) Mean number of ions (JVj) as a function of the position i in the channel, (b) Mean dimensionless electric potential 
{4>i) as a function of the position i in the channel. The channel is composed of L — 100 cells, each with iV_ = 25 fixed ions. 
The channel is submitted to the boundary conditions Nl = 45, Nn = 5, and 0l = 4>n- 




FIG. 2: (a) Spatial correlations Cjv(40,i) for the density fluctuations between site 40 and site i as a function of the position i 
in the channel, (b) Spatial correlations (7^(40, i) for the electric potential fluctuations between site 40 and site i as a function 
of the position i in the channel. The channel is composed of L — 100 cells, each containing iV_ = 25 fixed ions. The boundary 
conditions are JVl = 45, Nn = 5, and = <pR- 



which, by ergodicity, is equivalent to the ensemble average (A7j) = ^ A r i P(..., Ni, ...). We see that electroneutrality, 
(iVj) ~ N_, is well respected except for a small layer from the boundaries over a distance of the order of the screening 
distance A. The corresponding mean electric potential (</>,) is shown in Fig. [T]d. It is approximately linear throughout 
the channel, which corresponds to a constant electric field, except again for a small layer near the boundaries. 

Due to the electric interaction, the particle density fluctuations present spatial correlations along the channel. The 
correlations between cells are obtained by the time averages 

C N (i,j) ee (NiNj) - (N^ (N d ) - lim i f [N t (t) - (N,)} [N^t) - (Nj)] dt . (23) 

Fig. [2^, shows the effect of the electric interaction on the spatial correlations. We observe negative correlations between 
neighboring sites, showing that an excess charge at some location induces a repulsive effect on the neighboring sites. 
The negative correlations rapidly diminish with the separation between sites over a distance comparable to the Debye 
length. The presence of such correlations also shows that the particle distribution is not Poissonian even at equilibrium 
in contrast with population dynamics of chemical reactions network |23j . The electric potential correlation function 



C<f,(i,j) = ^>i<t>j) - (4>i) (<f>j) 



(24) 



G 



is depicted in Fig. |2j: 
and the boundaries. 



It presents long-range correlations with a linear decrease in both directions between the site 



IV. FLUCTUATIONS IN ION TRANSPORT 



In this section, we study the fluctuations in the transport process. For this purpose, we first introduce the fluctuating 
particle current j k counting the number of ions transferred between the cells k and k + 1. We denote by e k (s) = ±1 
the discrete charge transfer of an ion in the positive or negative direction between the cells k and k + 1 during the 
s th random transition occurring at the time t s [s k (s) = otherwise]. Accordingly, the fluctuating particle current is 
given by 



+oo 



Jk 



(t)= £ e k (s)6(t-t s ) 



(25) 



We notice that the currents ( 25 ) are so-called point processes composed of singular events occurring at random times. 



In the stationary state, the average particle current crossing the channel is time-independent and can be obtained 
from the stationary probability distribution as 



j=( Jk (t))= ]T pm,. 

N U ...,N L 



,N L )[W +k (N 1 ,...,N L )-W. k (N 1 ,...,N L )} 



(26) 



By current conservation, the average particle current is independent of the position in the channel, J = (jk) for all k. 

However, the experimentally measured electric current is composed of the particle and displacement currents ac- 
cording to Eq. ([3]). Indeed, even in a stationary state, fluctuations in the particle distribution at the mesoscopic level 
generate a fluctuating electric field that, in turn, generates a fluctuating contribution to the electric current. We thus 
calculate the change in the electric field associated with the random jumps of the charge carriers. The change in the 
electric field between the cells k and k + 1, E k+1 / 2 = — (^fe+i — &k)/Ax, is expressed as 



$i-#k+i + *k) =± 



"Ax 



k+lJ 



fc+l,Z+l 



k,l 



k,l+l 



(27) 



if the transition W±i is performed. The second equality is obtained by using Eq. (19) for the potential $'(..., Ni T 
1,-iVj+i ± 1,...) after the transition. The associated displacement current is given by ea times the change in the 
electric field, where a is the cross section of the channel (a Ax — AV). The displacement current associated with the 
transition W±i and measured between the cells k and k + 1 thus reads 




if k^l 
if k = l 



(28) 



where we used expression (A7) and the equality ea — aAx. Now, the contribution of the transition W±i to the total 
current is the sum of the particle current, ±e ii k = I and zero otherwise, and the displacement current (28) so that 
the location I contributes to the total current by the amount 



L + l 



m 



(29) 



as it should since the electric current must be divergence free at all times. Remarkably, all the transitions ±Z 
contribute equally to the electric current. The total current is finally obtained by summing the contributions from all 
the transitions ±Z: 



J(t) 



L 



t £*(*)• 



(30) 



1=0 



The sum runs over all the possible charge transfers since each charge displacement in the system modifies the electric 
field, inducing in turn a displacement current. Expression (30) can also be obtained from the Ramo-Shockley theorem 



11, 45 47 linking the current flowing in the external circuit to charge movement inside the system. By virtue of Eq. 
( 26 1 , the mean total current is related to the mean particle current according to 



(I) = eJ, 



(31) 
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FIG. 3: Power spectra of the particle current (evaluated at k = 5) and of the total current in units of e 2 (line with dots). The 
channel is composed of L — 20 cells and is submitted to the boundary conditions = 5, -/Vr, = 45, and (j>L = 0a- Each cell 
contains N- = 25 fixed ions. 



showing that the displacement current does not contribute on average, as it should. However, as shown below, the 
particle and total currents differ in the properties of their non-zero frequency fluctuations. 

Besides the mean currents, the noise characterization offers further information on the transport process [IT]. The 
power spectra of the currents are defined as the Fourier transform of the current correlation function: 



1 



+ 00 



W = ^ J e^[mm)-(0 \dt, (32) 

where £ is associated with either the particle current correlations, £ = j k , or to the electric current, £ = I. These 
power spectra can be obtained from an ensemble of long random trajectories according to 

S f ( U ) = T hm (||H| 2 ) - (0 2 8(w) (33) 

where £(ui) = L e~ %ut £(t)dt and where (•) denotes an average over several trajectories. For the particle current, we 
thus have 

jk{u) = I" e~ iut j k (t)dt = J2tk(s)e- iut » . (34) 



Similarly, the power spectrum of the total current ( 30 1 is obtained from 



!(<»)= I e "^w di = 7^iEE £ '( s ) e " ,wil ( 35 ) 



where the sums take into account every transition in the system, in accord with Eq. (30). 

In Fig. [3j we depict the power spectrum of the particle and total currents. We see that, under the same conditions, 
they are qualitatively different from each other, showing that the inclusion of the displacement current affects the 
fluctuations at positive frequencies and therefore must be properly taken into account. The power spectrum of 
the particle current increases with the frequency uj because the corresponding correlation function is negative. Such 
anticorrelations constitute another signature of the Coulomb repulsion between the ions, which reduces the probability 
of a definite transition event after jumping through the same section in the same direction. This frequency dependence 
of the power spectrum reflects the memory effects induced by the long-range Coulomb interaction on the random 
motion of the particles. On the other hand, the total current presents a strongly reduced noise at all frequencies 
because of the long-range correlations induced by the changes in the distribution of charges. We note that the power 
spectra reach a constant positive value in the high-frequency limit that results from the intrinsic randomness found 
at all times scales in such stochastic processes [35] . Precisely, in this limit, the spectra S k are given by the shot-noise 
formula 

S fc M= J2 Pst(N 1 ,...,N L ){W +k (N 1 ,...,N L ) + W- k (N 1 ,...,N L )}. (36) 

Nt,...,Nt 



In addition, in the high-frequency limit, the different noise processes become uncorrelated so that 

„2 l 



57(oo) 



(L + l) 



00 



(37) 



fc=0 



In the bulk of the channel away from the boundaries by a few Debye's lengths, the spectra Sk(uj) turn out to be 
uniform, i.e., approximately independent of the location k. This explains that, for a channel longer than Debye's 
length, we approximately find 



5*7(00) 



L + l 



S k (oo) . 



(38) 



as observed in Fig. [3] 

Furthermore, when the applied voltage V = $l — &r is sufficiently large we may neglect the backward transitions, 
W-k(Ni, Nl) ~ 0, in which case Eqs. (26) and (36) coincide so that 



5 fc (oo) = 



(I) 



(39) 



and 



5x(oo) = ^^(00) = ^ 



(40) 



by virtue of Eq. ( 37 ) . These expressions correspond to the Schottky noise formulae 



V. FLUCTUATION THEOREM FOR THE CURRENTS 



In this section, we focus on the fluctuations of the particle and total currents in the zero-frequency limit or, 
equivalently, in the long-time limit. More precisely, we will focus on the full probability distribution of the currents, 
which contain the information on the second-order noise properties as well as information on higher-order properties. 
In Ref. [34] it was shown in that the probability distribution of the particle current fluctuations may obey in the 
long-time limit, t 00, a fluctuation symmetry of the form 



jk(t')dt' 









a 










t Jo 



jk{t')dt' 



lAL 



(41) 



where A is the affinity (21 ). Here, this symmetry holds as a result of the thermodynamic properties of our stochastic 
description, detailed in Appendix|B| An illustration of this relation is given in Fig. [4^, where we depict the probability 
distribution function of the particle current, which is compared to the prediction of the fluctuation theorem (41 1. 



In addition, we here show that the fluctuation theorem also holds for the total current (30) exactly in the same 
form: 



Pi 



I(t')dt' = a 









a 










tJo 



I{t')dt' = -a 



^aAt 



(42) 



This symmetry can be rigorously shown by extending the demonstration of Ref. |34j to the present situation where 
any charge displacement in the system induces a non-local contribution to the total current, as shown by Eq. (30). 
The symmetry relation (42) is illustrated in Fig. [4Jd. The fluctuation relations (41) and (42) describe the large 
nonequilibrium current fluctuations and have important consequences on the linear and nonlinear response properties 

These results show that the fluctuation theorem is verified for systems of particles interacting with complex, long- 
ranged interaction, whereas most fluctuation relations were verified for locally interacting particles, for example with 
purely diffusive |5Ij or effusive |52j behavior, or with hard-core interaction [53] . For ion transport in channels, the 
fluctuation theorem was considered in the limiting case where the electric repulsion is strong enough to prevent the 
presence of more than one ion inside the channel |51[ I54H56| . 
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FIG. 4: (a) Probability distribution of the integrated particle current, G = J* jk(t')dt' , evaluated at k = 2. The comparison is 
performed with the fluctuation symmetry (41 1 for the negative part of the distribution (thick line with triangles), (b) Probability 
distribution of the integrated electric current, G = (1/e) J I(t')dt' . The comparison is performed with the fluctuation symmetry 
(42 l for the negative part of the distribution (thick line with triangles). The channel is composed of L — 10 cells, each containing 
N- — 25 fixed ions. The probability distributions are evaluated at time t = 273. The channel is submitted to the boundary 
conditions Ni, = -/Vr = 25 and (j>n — </>l = 0.05 so that the affinity takes the value A — 1/20. 



VI. CONCLUSIONS 



In this paper, we have introduced a mesoscopic description of ion transport in homogeneous channels on the ground 
of the laws of both electricity and statistical thermodynamics. The model is ruled by a master equation describing the 
spatial distribution of the discrete numbers of ions in the channel. Indeed, the stochastic aspects of the time evolution 
as well as the discrete nature of the ions are essential aspects of the behavior of matter at the mesoscopic scale. 
The model incorporates the self-consistent field generated by the fixed and mobile ions as described by the Poisson 
equation. Also, the model is shown to be consistent with thermodynamics and can be studied under equilibrium or 
nonequilibrium conditions. 

The stochastic description can be used to obtain the spatial correlations of particle density fluctuations in the chan- 
nel, which reveal the strong repulsion due to the electric interaction for fluctuations that depart from electroneutrality. 
Moreover, we observe long-range correlations in the electric potential fluctuations. 

Another quantity of fundamental interest in the study of ion transport in channels is the experimentally measured 
total current, which is composed of the particle current plus the displacement current arising from the temporal 
variations of the electric field. At the macroscopic level, the displacement current vanishes in a stationary state. On 
the other hand, at the mesoscopic level, the molecular fluctuations induce long-ranged electric field fluctuations which 
result in a reduced noise spectrum at all positive frequencies. The particle and total currents thus differ in their 
fluctuation properties and the present approach allows us to assess the importance of such effects. Furthermore, we 
have shown that the fluctuation theorem for the currents holds in this more general situation as well, that is when 
the displacement current is taken into account in the total current. 

In summary, we have shown how to describe ion transport in homogeneous channels at the mesoscopic level. The 
extension of the present considerations to the transport of several ionic species in one channel and to the case of 
heterogeneous channels of relevance to biological systems is possible on the basis of the present results. 

Acknowledgments. D. A. is grateful to the F. R. S.-FNRS for financial support. This research is financially 
supported by the Belgian Federal Government (IAP project "NOSY") and the "Communaute frangaise de Belgique" 
(contract "Actions de Recherche Concertees" No. 04/09-312). 



Appendix A: Solution of the channel electrostatics 



The electric potential $j is defined on each cell i = 1,...,L and obeys the discretized Poisson equation (11) (above 
written in terms of the dimensionless potential <f> — —e^/k^T). This linear system must be solved at each time the 
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particle distribution changes and must satisfy the boundary conditions 

$ = $ L , $ L+1 =$ R . 
This linear system of equations can be written in matrix form as 

C $ = Z 



with the vectors 



and 



$ T = ($ 1; ...,$ L ) 



Z T = e(JVi-JV_,...,iV i -JV_) + a(*L,0,...,0 ) $ R ) 
= Q + a($L,0,...,0,$R), 



(Al) 
(A2) 
(A3) 

(A4) 



where T denotes the transpose operation and where we introduced the charge vector Q of components Qk — e(Nk — 
N- ) . Note that the boundary conditions appear on the first and last components of the vector Z multiplied by the 
quantity a = (eA r _)/(n_ Ax 2 ), which has the units of a capacitance. The L x L symmetric matrix C reads 



a 



( 2 -1 
-1 2 -1 

-1 2 

-1 



-1 

2 -1 
-1 2 J 



(A5) 



It is invcrtible so that we have 



with 



$ = C 1 Z 



-^(L + l-j) if i<j 
-^(L + l-i) if i>j. 



(A6) 



(A7) 



The matrix C 1 is symmetric, as it should. 

The electrostatic energy U stored in the system is given by 



U= I$T. C . $= I Z T. C -1. Z 

2 2 



(A8) 



The change in electrostatic energy associated with the transition of an ion from cell i to cell i + 1 is thus given by 

AU i>i+1 - l - (Z' T • CT 1 • Z' - Z T • C- 1 • Z) (A9) 

where Z' k = Zk — eSks + e$k,i+\ characterize the change in the charge distribution resulting from the transition 
i-^i + 1. Developing this expression yields 



2 

AC4 m = e(<i> m - $0 + 6 - f (C- 1 ) -2(C- X ) +(C" 1 ) 

Z L \ / i,i V / i.i+1 \ /»+!,( 



+ 1 



(A10) 



where we used the symmetry of the coefficients (C 1 ) = (C 1 ) as well as the expression $j ■ = J2k (C 1 ) ^fe 

for the electric potential in cell j. We see that the change in electrostatic energy Ai7j.,- + i only depends on the initial 
voltage difference plus a term independent of the charge state of the system. In the case where the transition involves 
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one of the reservoirs, say the left one, we have $l = $l and the change in electrostatic energy for a transition from 
the left reservoir to cell 1 reads 



AC/ 0jl = e ($ 1 -$ L ) + ^(c- 1 ) i 



and similarly for transitions involving the right reservoir. By virtue of expression (A6), we see that 

*i = E ( c_1 ) , z 'k = E ( c_1 ) , ^ - e **.< + e5k ^ = *i - e ( c_1 ) + e ( c_1 ) 



(All) 



(A12) 



which is Eq. (191 for the transition i — > i + 1. Accordingly, the energy differences ( |A10 1 and (All) can be expressed 
in terms of the initial and final voltage differences Vi and V[ between the cells i and i + 1 so that we recover Eq. ( 16 ) 
for the electrostatic energy difference. Note that Eq. ( 16 ) also holds for inhomogeneous channels with e = e(x) and 
ri_ = ri-(x). 



Appendix B: Macroscopic affinities of the stochastic model 




FIG. 5: Graph associated with the random process of the master equation ( 12 I over a chain of length L — 2. The depicted 
cycle corresponds to the transfer of an ion from the left to the right reservoir. The orientations of the edges are chosen so that 
a movement from left to right is counted positively. 



In this appendix, we detail the thermodynamic properties of the stochastic model introduced in Sec 
a system ruled by the master equation (12 1, a graph G is associated as follows [ 



each state 



TIT] For 

of the system 



corresponds to a vertex or node while the edges represent the different transitions lu ^ lj' allowed between the states. 
The stochastic ion transport of Sec. Ill is defined in terms of the populations (N%, ...,Nl) of mobile ions inside the 



channel. Accordingly, each such configuration defines a node of the graph and each transition ±i corresponds to one 
edge. As an example the graph associated with the case L — 2 is drawn in Fig. [5] 

As a matter of fact, cyclic trajectories play a special role as they link the macroscopic thermodynamic properties 
to the transition rates of the mesoscopic description [22 [21]. In the case L = 2 corresponding to Fig. [5j we see that 
any cyclic path in the graph can be decomposed as a linear combination of two types of cyclic trajectories: 



Cl 



(Ni,N 2 ) -> (Ni + 1, N 2 ) -> (Ni + 1, N 2 + 1) -> (N U N 2 + 1) -> (N u N 2 ) 



(Bl) 



and 



c 2 = (N lt N 2 ) -> (JV X + 1, N 2 ) ~> (N lt N 2 + 1) -> (N 1: N 2 ) 



(B2) 
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As shown by Schnakenberg [55], the macroscopic affinities associated with these cyclic paths can be obtained from 
the transition rates by calculating the quantities 



A(ci) = In 



W +0 {N 1 ,N 2 )W- 2 {N 1 + 1, jV 2 )W_ (JV 1 + 1, JV 2 + l)W +2 (iVi, jV 2 + 1) 
W_ (iVi + l,^V 2 )iy +2 (7V 1 + 1, AT 2 + i)W +0 (JVi,JV a + l)W_a(JVi,JVa) 



and 



A(c 2 ) = In 



W +0 (JVi, JV 2 )W+i(JVi + 1, JV 2 )W +2 (JVi, JV 2 + 1) 

w_ (JVi + i,jyr 2 )w_i(JVi,j\r2 + i)w_ 2 (jvi,jv 2 ) 



(B3) 



(B4) 



These expressions are obtained as the ratio between the product of the transition rates along the cyclic path in one 
direction over the product of the transition rates along the cyclic path in the reversed direction. Calculating the 
affinities (B3) and (B4| with the transition rates (14) yields 



A(a) = and A(c 2 ) = In 



Nq 



L + l 



(B5) 



The first one always vanishes whereas the second is the macroscopic affinity (21 ), irrespectively of the initial configu- 
ration (Aq, Nl). Note that A(c 2 ) only involves externally controlled parameters, as it should. This results from the 
fact that only the second cycle involves the transport of an ion across the channel. The extension to larger channels, 
L > 2, is straightforward. 

According to a result of Kolmogorov |57j , a Markov process is at equilibrium if and only if the affinities along the 
cycles vanish. As shown by Eqs. (B5) this occurs when condition fcty is fulfilled. 
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